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We investigate and discuss how hydrogen behaves at the edges of a graphite sheet, in par- 
ticular the armchair edge. Our density functional theory-based calculations results show that, 
in contrast to the zigzag edge [cf., e-J. Surf. Sci. Nanotech. 2 (2004) 77], regardless of orienta- 
tion, there is an activation barrier hindering H2 dissociation at the armchair edges. And once 
they do get dissociatively adsorbed at the armchair edges, we find that it would be extremely 
hard to desorb the H from their adsorption sites at the armchair edges. Furthermore, we also 
found that, consistent with our earlier conclusions [cf., J. Phys. Soc. Jpn. 72 (2003) 1867], it 
is unlikely that we would find a whole H2 in between plain graphite sheets. 
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Carbon-based nanomaterials, e.g., carbon nanotubes 
(CNTs) and graphite nanofibers (GNFs), have been at- 
tracting much attention because of their purported po- 
tentials as materials for gas storage. Earlier reports of 
high hydrogen uptake of these materials make them at- 
tractive as hydrogen storage devices in fuel-cell-powered 
electric vehicles.^' ^ However, although various theoreti- 
cal and experimental studies have been made (cf., 3-18, 
and references therein) , and initial reports seem to show 
promising and spectacular results (more than 60 wt% in 
some cases, which is way above the minimum 6.5 wt% 
hydrogen storage density requirement for economically 
feasible mobile application of hydrogen fuel at room tem- 
perature), research in this field remains contentious, as 
later studies reveal trends significantly less encouraging 
than were first contemplated.^^' This situation is in- 
dicative of the necessity for more systematic investiga- 
tions. On top of that, there is a necessity to develop 
a microscopic picture of the mechanism underlying hy- 
drogen adsorption onto, absorption into, and desorption 
from carbon related materials. With this in mind, we 
have done, and are still continuing, a series of studies 
on the interaction of hydrogen with various surfaces of 
graphite. ^^"^^ 

Invoking the density functional theory (DFT), we cal- 
culated the potential energy surfaces (PES) relevant to 
the dissociative adsorption of H2 onto a graphite sheet. 
Our results show that the reconstruction of the C atoms 
plays an important role in understanding the H2-graphite 
surface interaction. We observed a lowering of the disso- 
ciation barrier by ca. 1 eV, as a result of the relaxation 
outwards of the C atoms (changing its structure from 
sp^ to sp'^-like), to meet the incoming H2. However, this 
is not enough to allow H2 to dissociate, as the dissocia- 
tion barriers still remains high at ca. 3 eV. Furthermore, 
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because the C atoms have to move out and meet the 
incoming H (H2), this becomes difficult for H2 dissocia- 
tion, esp., when the C atoms involved are very close to 
each other (ca. 1.42 A). The PES for H2 dissociation on 
the surface plane of a graphite sheet is thus azimuthally 
corrugated. The relaxation of the C atoms also stabilizes 
the H-graphite surface interaction, with a corresponding 
binding energy of 0.67 eV. We were later informed^^ that 
the theoretical results quantitatively agree with recent 
experimental results. 

These results also give us an idea as to how hydro- 
gen would behave when it interacts with a single-walled 
carbon nanotube (SWCNT).^^ As expected from earlier 
results, ^^'^'^ a SWCNT with a small diameter adsorbs 
H quite easily. However, as the diameter increases, the 
potential energy features approach that for hydrogen in- 
teraction with the surface plane of a graphite sheet. It 
becomes more and more difficult to dissociate H2 as the 
SWCNT diameter increases. 

We also considered the possibility of finding hydrogen 
in between graphite layers. Again, the reconstruction 
of the C atoms played an important role in determining 
what stable configurations hydrogen will assume once 
found inside/between graphite layers — the H atom as- 
sumes a position ca. 1 A from one C atom on one sheet 
and ca. 2 A from the C atom on the other sheet. On the 
other hand, it would be difficult to dissociate H2 in be- 
tween graphite layers, let alone find them as H2 between 
the graphite layers. Again, we found amazing agreement 
between our results and recent observations from neu- 
tron diffractiou"'^^ and thermal desorption^* studies on 
deuterated nanographite. 

Having found that it would be difficult for the H atom 
to reach inside/between the graphite layers through the 
hexagonal holes at the surface planes, ^^'^'^ and that it 
would also be difficult to find H2 inside/between graphite 
layers,^'* we considered the next probable surfaces on 
graphite, viz., the edges. Our calculation results show 
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Fig. 1. (Upper Panel) H2 configuration with respect to the arm- 
chair edge, the H-H bond is oriented parallel to and immediately 
above one of the graphite sheets as indicated. H2 is constrained 
to dissociate over the indicated site in a parallel orientation and 
a planar geometry. (Lower Left Panel) Contour plot of the 
calculated potential energy surface (PES) for H2 dissociation at 
the indicated site on the armchair edge of graphite, as a function 
of the H-H separation r, and the normal distance Z of the H2 
center-of-mass from the armchair edge. The dashed-line traces 
the 'path of least resistance/potential (reaction path). The in- 
terlayer distance between graphite sheets is ca. 3.5 A. (Lower 
Right Panel) Calculated potential energy along the reaction 
path, i.e., along dashed-line indicated on the PES shown on the 
Lower Left Panel. Energies are given in electronvolts rela- 
tive to the values at {Z = 00, r 0.74 A). Contour spacing is 
Si 0.2 eV. 



that the zigzag edge is very reactive, and it is possible 
to dissociatively adsorb H2 at the zigzag edge without 
any activation barrier hindering the reaction. We con- 
tinue our study here, and consider the interaction of hy- 
drogen with the other remaining edge of graphite, viz., 
the armchair edge. 

We perform DFT-based total energy calculations, 
using the plane waves and pseudopotentials.^* The 
graphite supercell used in the calculations consists of 
two graphite sheets/layers, which are separated by ca. 
3.5 A and shifted relative to each other as shown on the 
Upper Panels of Figs. 1 to 5. Each graphite layer consists 
of 21 carbon atoms with a nearest neighbor distance of 
1.42 A . We have chosen a sufficiently large (graphite) 
supercell to avoid interactions between the H atoms (H2 
molecules) in the neighboring supercells. The electron- 
ion interaction is described by optimized ultrasoft pseu- 
dopotentials,^^ and the Kohn-Sham equations are solved 
using plane waves with kinetic energies up to « 476 eV. 
The surface Brillouin zone integration is performed using 
the special-point sampling technique of Monkhorst and 
Pack (with 4x4x1 sampling meshes). For the ex- 
change correlation energy, we adopt the generalized gra- 
dient approximation (GGA).'^-'^''^'^ No significant change 
in the numerical results was observed when we increased 
the kinetic energy cutoff and the number of sampling 
points. 

Figures 1 and 2 (Lower Left Panels) show the calcu- 
lated potential energy contours for H2 dissociation on 
the armchair edge of graphite, as a function of the H- 
H separation r, and the normal distance Z of the H2 




Fig. 2. (Upper Panel) H2 configuration with respect to the arm- 
chair edge, the H-H bond is oriented parallel to and immediately 
above one of the graphite sheets as indicated. H2 is constrained 
to dissociate over the indicated site in a parallel orientation and 
a planar geometry. (Lower Left Panel) Contour plot of the 
calculated potential energy surface (PES) for H2 dissociation at 
the indicated site on the armchair edge of graphite, as a function 
of the H-H separation r, and the normal distance Z of the H2 
center-of-mass from the armchair edge. The dashed-line traces 
the path of least resistance/potential (reaction path). The in- 
terlayer distance between graphite sheets is ca. 3.5 A. (Lower 
Right Panel) Calculated potential energy along the reaction 
path, i.e., along dashed-line indicated on the PES shown on the 
Lower Left Panel. Energies are given in electronvolts rela- 
tive to the values at {Z = 00, r Ri 0.74 A). Contour spacing is 
0.2 eV. 



center-of-mass from the armchair edge. The H2 is ini- 
tially oriented with the H-H bond parallel to one of the 
graphite sheets, and the dissociation is constrained in a 
parallel orientation with respect to the armchair edge of 
one of the graphite sheets as indicated in the Upper Pan- 
els of Figs. 1 and 2. H2 dissociation is constrained in a 
parallel orientation with respect to the armchair edge, 
and a planar geometry. Energies are given in electron- 
volts relative to the values for the case when H2 is in the 
gas phase, far {Z = 00) from the graphite. The interlayer 
distance between graphite sheets is ca. 3.5 A. 

We can see that there is an activation barrier hindering 
H2 dissociation on the armchair edge of graphite. There 
is ca. a 1.7 eV barrier hindering H2 dissociation in the 
configuration depicted in Fig. 1 (cf., corresponding Lower 
Right Panel), and ca. 0.6 eV hindering H2 dissociation 
for the configuration depicted in Fig. 2 (cf., correspond- 
ing Lower Left Panel). Furthermore, assuming that we 
do succeed in dissociating H2 at the indicated sites, we 
encounter the problem of how to desorb these H atoms, 
as they stick strongly to the corresponding edges, with 
an adsorption energy ranging from ca. 4 to 5 eV. 

Reorienting the H2 in the configuration depicted in 
Figs. 3 and 4 (Upper Right Panels), such that the H- 
H bond is now perpendicular to each of the graphite 
sheets, did little to improve the situation. For the con- 
figurations depicted in Fig. 3, the activation barrier still 
remains high at ca. 0.5 eV (Fig. 3, Lower Right Panel). 
And, although we observe a significant lowering of the 
activation barrier to a more or less manageable ca. 0.3 
eV in the configuration depicted in Fig. 4 (Lower Right 



H2 Dissociative Adsorption at the Armchair Edges of Graphite 



3 




oO-O*— »0-0«— <K>"0®— f 
(I <■ H: 

0— wO"0* — !>0-0*' — »o-o« 




Fig. 3. (Upper Right Panel) H2 configuration with respect to 
the armchair edge, the H-H bond is oriented perpendicular to 
and in between the two graphite sheets as indicated. H2 is con- 
strained to dissociate over the indicated site in a parallel orien- 
tation and a planar geometry. (Left Panel) Contour plot of the 
calculated potential energy surface (PES) for H2 dissociation at 
the indicated site on the armchair edge of graphite, as a function 
of the H-H separation r, and the normal distance Z of the H2 
center-of-mass from the armchair edge. The dashed-line traces 
the path of least resistance/potential (reaction path). The in- 
terlayer distance between graphite sheets is ca. 3.5 A. (Lower 
Right Panel) Calculated potential energy along the reaction 
path, i.e., along dashed-line indicated on the PES shown on the 
Left Panel. Energies are given in electronvolts relative to the 
values at (Z = 00, r Ri 0.74 A). Contour spacing is 0.2 eV. 
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Fig. 4. (Upper Right Panel) H2 configuration with respect to 
the armchair edge, the H-H bond is oriented perpendicular to 
and in between the two graphite sheets as indicated. H2 is con- 
strained to dissociate over the indicated site in a parallel orien- 
tation and a planar geometry. (Left Panel) Contour plot of the 
calculated potential energy surface (PES) for H2 dissociation at 
the indicated site on the armchair edge of graphite, as a function 
of the H-H separation r, and the normal distance Z of the H2 
center-of-mass from the armchair edge. The dashed-line traces 
the path of least resistance/potential (reaction path). The in- 
terlayer distance between graphite sheets is ca. 3.5 A. (Lower 
Right Panel) Calculated potential energy along the reaction 
path, i.e., along dashed-line indicated on the PES shown on the 
Left Panel. Energies are given in electronvolts relative to the 
values at (Z = 00, r 0.74 A). Contour spacing is fa 0.2 eV. 



Panel) , there still remains the problem of how to desorb 
the H atoms. The adsorption energy still remains high 
at ca. 4 to 4.5 eV (Figs. 3 and 4, Lower Right Panels). 

We also checked to see whether or not it would be 
possible to insert the H2 into/between the graphite lay- 
ers, as depicted in Fig. 5 (Upper Right Panel), and then 
dissociate them. Consistent with our earlier results,^"* 
we find that it would be difficult to dissociate H2 in be- 
tween graphite layers, let alone find them as H2 between 
the graphite layers. 

These results may be compared to the activated dis- 
sociative adsorption on a graphite sheet^^'^^'^^ and non- 
activated dissociative adsorption on the zigzag edge.^^ 
The spectacular change in reactivity mainly lies in the 
difference in structure of the C atoms at these surfaces. 
On a graphite sheet, all of the C atoms are bonded to 
three neighboring C atoms, with sp^-structure. The ro- 
bustness of the sp^-structure prevents the C atoms on 
the graphite sheet from forming new bonds. Recall that 
a C atom must first form an sp^-like structure before 
it can accomodate a H atom, directly impinging on the 
graphite sheet, or one coming as a product of H2 dis- 
sociation.^^' ^■^■^'^ Compared to the graphite sheet, the 
C atoms at the outermost edges (armchair and zigzag) 
are bonded to only two neighboring C atoms. This makes 
the edges more reactive (more susceptible to forming new 
bonds) to H adsorption and/or H2 dissociation, as com- 
pared to the graphite sheet. The difference in reactivity 
between the armchair and the zigzag edge lies in the way 
the C atoms at the edges are bonded. As mentioned ear- 
lier, each of the outermost C atoms at the armchair edge 



has two neighboring C atoms, one of which is also located 
at the armchair edge. In effect, every reaction occurring 
at an armchair edge, e.g., H adsorption or H2 dissocia- 
tion, involves an armchair edge C-C pair. When one of 
the C atoms in the armchair edge C-C pair forms a new 
bond with the incoming H atom, the other remaining C 
atom (of the C-C pair) must somehow compensate for 
the resulting change in covalency. Although a C atom 
located at the outermost part of a zigzag edge is also 
bonded to two neighboring C atoms, neither of these 
neighboring C atoms are located at the edge. The C 
atoms at the zigzag edge are, thus, already in a favor- 
able configuration, and ready to welcome the incoming 
hydrogen. (Here, we restrict the discussion to geomet- 
ric effects and defer the discussion on the corresponding 
electronic effects. '^^) 

In summary, based on the density functional theory, 
we investigate and discuss how hydrogen behaves at the 
edges of a graphite sheet, in particular the armchair edge. 
We found that, regardless of site, there is always an ac- 
tivation barrier hindering II2 dissociative adsorption — 
a spectacular difference compared to the very reactive 
zigzag edge.^^ Furthermore, we also found that it was 
difficult to insert the whole H2 through the armchair 
edge, in agreement with our earlier findings in our stud- 
ies of II2 inside/between graphite layers.^"* Based on the 
results we have so far gathered, we conclude that to be 
able to use carbon nanomaterials as a means to store 
hydrogen, the crucial steps would be to dissociate hy- 
drogen first, and then (somehow) induce them to stick 
to the carbons on each sheets. The results we present 
suggest the possible utility of the zigzag edge over the 
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Fig. 5. (Upper Right Panel) H2 configuration with respect to 
the armchair edge, the H-H bond is oriented parallel to and in 
between the two graphite sheets as indicated. H2 is constrained 
to dissociate over the indicated site in a parallel orientation and 
a planar geometry. (Left Panel) Contour plot of the calcu- 
lated potential energy surface (PES) for H2 dissociation at the 
indicated site on the armchair edge of graphite, as a function 
of the H-H separation r, and the normal distance Z of the H2 
center-of-mass from the armchair edge. The dashed-line traces 
the path of least resistance/potential (reaction path). The in- 
terlayer distance between graphite sheets is ca. 3.5 A. (Lower 
Right Panel) Calculated potential energy along the reaction 
path, i.e., along dashed-line indicated on the PES shown on the 
Left Panel. Energies are given in electronvolts relative to the 
values at (Z = 00, r fs! 0.74 A). Contour spacing is w 0.2 eV. 



armchair edge as a reaction channel to reaUze efficient 
efficient hydrogen-storing carbon nanomaterials. 
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